Ordering and Fluctuation of Orbital and Lattice Distortion 
in Perovskite Manganese Oxides 
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Roles of orbital and lattice degrees of freedom in strongly correlated systems are investigated to 
understand electronic properties of perovskite Mn oxides such as Lai-ajSrxMnOa. An extended 
double-exchange model containing Coulomb interaction, doubly degenerate orbitals and Jahn- Teller 
coupling is derived under full polarization of spins with two-dimensional anisotropy. Quantum fluc- 
tuation effects of Coulomb interaction and orbital degrees of freedom are investigated by using the 
quantum Monte Carlo method. In undoped states, it is crucial to consider both the Coulomb inter- 
action and the Jahn- Teller coupling in reproducing characteristic hierarchy of energy scales among 
charge, orbital-lattice and spin degrees of freedom in experiments. Our numerical results quantita- 
tively reproduce the charge gap amplitude as well as the stabilization energy and the amplitude of 
the cooperative Jahn- Teller distortion in undoped compounds. Upon doping of carriers, in the ab- 
sence of the Jahn- Teller distortion, critical enhancement of both charge compressibility and orbital 
correlation length is found with decreasing doping concentration. These are discussed as origins of 
strong incoherence in charge dynamics. With the Jahn- Teller coupling in the doped region, collapse 
of the Jahn- Teller distortion and instability to phase separation are obtained and favorably com- 
pared with experiments. These provide a possible way to understand the complicated properties of 
lightly doped manganites. 



PACS numbers; 71.30.-|-h, 71.27.+a, 71. 20. Be, 71.38. -l-i 



I. INTRODUCTION 
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Perovskite Mn oxides show a variety of interesting 
physics depending on chemical composition and tem- 
perature. Pioneering work in 1950's revealed fascinat- 
ing phase transitions froin_ap. antiferromagnetic insulator 
to a ferromagnetic metaLBiZl The transitions take place 
through a change in the chemical composition which 
effectively introduces mobile holes in the systems. In 
1970, a large negative magnetoresistance was discovered 
near the transition from the ferromagnetic metai-to a 
paramagnetic state with increasing temperaturcHiJ Re- 
cently, this so-called colossal magnetoresistance has at- 
tracted renewed interest not only for understanding ba- 
sic mechanispas but also for a possibility of technological 
applications .E3 

The basic physics in these systems has been exam- 
ined by the so-called double-exchange mechanism. El In 
undoped states, each Mn'^+ ion has four electrons in 
the {t2gY{eg)^ configuration. Doped mobile holes are 
introduced to the twofold Cg orbitals. These situations 
have been examined based on the j-simplest version of 
the double-exchange (DE) model JjQlZlOTi3 The model 
consists of noninteracting conduction electrons in a non- 
degenerate orbital and localized spins ferromagnetically 
coupled to this conduction band through the Hund's rule 
couplingE The explicit form of the Hamiltonian is given 
by 



where c|^(cicr) creates (annihilates) a cr-spin electron at 
site i; CTj is the spin operator for conduction electrons 
whereas Si denotes the localized spin at site i. This 
model successfully explains the metallic ferromagnetism 
under hole doping, the transition to paramagnetic phase 
with increasing temperature and the colossal magnetore- 
sistance. 

Recently, intensive experimental studies have revealed 
an importance of interplay among charge, spin_orbital 
and lattice degrees of freedom in the Mn oxides £j These 
systems show a Mott insulating or a charge-ordered state 
concomitant with spin, orbital and lattice orderings; a 
metallic state with different spin and orbital orderings; 
and various transitions between them. These novel phe- 
nomena strongly indicate a necessity for studies beyond 
the 'simple' DE model given by Eq. (|l|). 

In particular, we discuss here the following experi- 
mental results which are difficult to understand within 
the 'simple' DE model: In undoped states, these ma- 
terials are insulating with a charge gap of the order of 
leV. At low temperatures, this state shows spin, or- 
bital and lattice orderings with strong two-dimensional 
(2D) anisotropy. The magnetic structure is called A-type 
and consists of ferromagnetic planes stacked antiferro- 
magnetically. For instance, in LaMnOa, the charge gap 
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amplitude Ac is estimated around leV, the tran- 
sition to the 2D orbital-ordered state with lattice dis- 
tortion (cooperatiifft-Jahn- Teller (JT) distortion) occurs 
at Tc ~ 780Kjj|'ElErEl and thcpjd^ype antiferromagnetic 
phase appears at Tn — 140K.ala't3 These indicate a hier- 
archy structure of energy scales: Ac ^ Tc ^ 7n- 

When we dope holes, the materials show com- 
plex transitions at low temperatures; collapse of 
the JT diS|tji|ff|t|bH^_magnetic and metal-insulator (MI) 
transitions-Era'clEllljCj lJ Doped holes generally induce 
a spin canting; the spin direction in the 2D ferro- 
magnetic planes tilts and a net ferromagnetic mo- 
ment is induepd--, without destroying the cooperative 
JT distortion.cHalla'cSEa Recently, neutron scattering ex- 
periment reported a charge or polaron ordering with 
2D anisotropy around S = 0.125 in the case of 
Lai_2;Sr2;Mn03 (Sr co n c en tration x corresponds to the 
hole concentration (5).E^ED In this phase, the coopera- 
tive JT distortion is not observed whereas the spins are 
canted. The distortion collapses at 5 ~ 0.1. The system 
remains insulating through the spin-canted state and fi- 
nally at 5 ~ 0.16, becomes the ferromagnetic metal. 

An important and puzzling feature appears in the fer- 
romagnetic state near the MI transition; in the case of 
Lai-a^Srj^MnOs, for S ^ 0.16. The metal exhibits strong 
charge incoherence at low temperatures where spin de- 
grees of freedom are frozen up due to the perfect po- 
larization. A strong mass enhancement is suggeste d hy . 
a small Drude weight D in the optical conductivityn3'E2l 
and a small discontiauity at the Fermi edge in the pho- 
toemission spectraJ23 while there is no conacicuous en- 
hancement in the specific heat coefficient 7.EJ 

In order to consider the above experimental aspects, 
several approaches beyond the 'simple' DE model (|l|) 
have been explored. One is mean-field approximation. 
The undoped states have been investigated by effective 
Hamiltonian derijvsd-irom a strong coupling limit of elec- 
tron correlationsj23"E3 by extended_DE models with both 
spin and orbital degtees of freedom,E3 and by models con- 



taining oxygens. I 



I These mean-field results commonly 



suggest energetic degeneracy between different spin and 
orbital structures. The importance of JT coupling or 
i2g-spin interaction has been discussed to lift this degen- 
eracy and to reproduce the realistic state. The following 
points, however, are missing in these calculations: They 
have not seriously considered the experimental energetic 
hierarchy, Ac > Tc > Tn- Although the JT distortion 
appears in much higher energy scale than the spin order- 
ing temperatures, the optimization of the ground state 
including phonons has not been successfully done. The 
model parameters, for instance, the interaction strength 
or the JT coupling they have chosen, are rather scattered 
and the consensus of parameter values is not reached. 
Among others, it should be noted that, as suggested in 
some degeneracy of various phases in these calculations, 
the magnetic and orbital structure must be under strong 
quantum fluctuations. 

There is another approach, first-principles calculation. 



which is applied mainly to undoped states. The local- 
density approximation (LDA) and the local-spin-density 
approximation (LSDA) conclude a gapless metal for the 
undoped state when the JT distortions are neglected. 
On the assumption of the experimental lattice structure, 
LSDA gives the antiferromagnetic insulator, although 
the charge gap amplitujdepis insufficient in reproducing 
the experimental result .LJO These are a consequence of 
the general tendency to underestimate electron correla- 
tion effects in these approximations. LDA-I-U calcula- 
tions, in which the effects of strong correlation are taken 
into account by a Hartree-Fock scheme, have also been 
done under the experimental lattice structure and repro- 
duced realistic values Afi-the charge gap amplitude and 
the magnetic moment.ocll The improvement reached in 
LDA-I-U with the JT distortion implies that not only the 
JT coupling but also the strong correlation play an im- 
portant role in these materials. The correct description 
of the experimentally-observed ground state, especially 
the large charge-gap amplitude, has not been reached ab- 
initio. We note that it is rather hard to examine quantum 
fluctuation effects by the available band-structure calcu- 
lations. 

In contrast to the undoped states, doping effects have 
been less investigated theoretically. Especially, the JT 
distortions in doped regions have not been seriously ex- 
amined. A systematic study on the doping effects is de- 
sired to discuss the collapse of the distortions and the 
charge or polaron ordering in experiments. 

In the ferromagnetic metal, the strong charge incoher- 
ence has been discussed beyond the 'simple' DE models. 
Since the spins are perfectly polarized at low tempera- 
tures, this incoherence strongly suggests an importance 
of neglected elements in the model. Many candidates of 
the drivingjmechanism have been discussed such £ys-|JHi 
fluctuation,E3 interorbital hoppiiig of the eg orbitals,E2rc3 
and strong Coulomb interaction.L3 Because the charge in- 
coherence conspicuously appears near the MI transition, 
it is important to fully understand the undoped insulat- 
ing states and to consider doping effects on them in a 
unified framework for treating quantum fluctuations. 

In the present paper, ground state properties of the 
extended DE model containing the Coulomb interaction, 
twofold degeneracy of Cg orbitals, and the JT coupling are 
discussed. The undoped and doped states as well as criti- 
cal properties of the MI transition are examined in a uni- 
fied framework. A part of thexesults has been published 
in the previous publication. cJa In this paper, more sys- 
tematic analyses are presented including the JT effects 
in hole-doped regions. As mentioned above, the real ma- 
terials show strong 2D anisotropy at and around the un- 
doped state, which suggests the importance of quantum 
fluctuations. The JT coupling to some extent suppresses 
these quantum fluctuations and plays a crucial role to re- 
produce experiments. In this work, the quantum fluctu- 
ation effects are investigated in the 2D anisotropic struc- 
ture beyond the mean- field results; we use here an unbi- 
ased approach by employing the quantum Monte Carlo 
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method.caEll Effects of tlie JT distortion are also studied 
by optimizing the ground state of the Hamiltonian with 
an electron-phonon couphng. Physical properties in the 
optimized ground state are calculated and compared with 
the experimental results with emphasis on roles of the hi- 
erarchy structure in energy scales. We determine values 
of model parameters appropriate for the Mn compounds 
and compare with those in other theoretical calculations. 
For doping of carriers, critical enhancement of both or- 
bital correlation length and charge compressibility is ob- 
tained in the absence of the JT distortion. The results 
are examined to understand the puzzling feature of the 
incoherent metal in experiments. We also show that the 
JT distortion rapidly becomes unstable with increasing 
doping. Phase separation takes place in our model when 
we take account of the JT coupling in the doped cases, 
which implies an instability to charge ordering in real 
materials. 

This paper is organized as follows: In Sec. ||, we de- 
scribe derivation of our model and discuss the parameters 
used in the following discussions. Numerical results are 
presented in Sec. III. The numerical method and con- 
ditions are also described. All the results are discussed 
in comparison with experiments as well as with other 
theoretical studies in Sec. IV. Section ^ is devoted to 
summary. 



where Hhop = J2ij.,u,^',at'ij' 4i^a<^j'^''^ which describes 
electron hoppings; Hint denotes on-site electronic inter- 
actions within doubly-degenerate Cg orbitals; Hk is sim- 
ilar to the second term in Eq. ([l]) which describes the 
Hund's rule coupling between electrons in the Cg orbitals 
and localized spins in the t2g orbitals; and Ti.^ controls 
total electron density by the chemical potential fi. The 
indices i,j = 1, • • •,iVs and ^,1/' = 1,2 denote the lat- 
tice sites of Mn atoms and the orbitals, respectively. We 
take ly = 1(2) for the dj.2_y2{d:}z'^_r2) orbital. Although 
the t2g-spin interaction such as rS,ir^j has been consid- 
ered in several previous studies,c3i23 it is ignored in Eq. 
(^) because the e nerg y scale is small. This point will be 
discussed in Sec. |II D| . 

Following several previous studies ,eH3 we consider the 
limit of the strong Hund's-rule coupling, K\S^'^^\ :s> t'^-' , 
which may be appropriate for situations in many per- 
ovskite manganites. Since (Ti aligns parallel to Si in each 
site in this limit, the hopping integral between sites i and 
j is renormalized by a relative angle of the spins at sites 
i and j. Then the model (H) is simplified into a spinless 
fermion model as 



(4) 



II. MODEL 



In order to examine the above-mentioned experimental 
results beyond the 'simple' DE model defined in Eq. (||), 
we take account of on-site Coulomb interaction, twofold 
degeneracy of conduction bands, anjd_JT coupling. Our 
Hamiltonian consists of three terms O 



(2) 



The first electronic term Tid is derived from the extended 
DE model which contains the Coulomb interaction and 
two orbi tal d egrees of freedom of eg electrons as described 
in Sec. II A in detail. The second and third terms de- 



note the electron-phonon interaction (JT coupling) and 
the elastic energy of phonons, respectively. Phonons are 
treated in a classical way as mentioned in Sec. |II B| . Sec- 
tion fl C| describes details of the mo del p arameters used 

the derivation 



II D 



in the following calculations. In Sec. 
of our model is examined to clarify energy scales and 
doping regions where our model can be applied. 



A. Electronic part 



The extended DE model containing twofold degeneracy 
F Cg orb 
given by 



of eg orbitals as well as on-site Coulomb interaction is 



^DE ~ ^hop 



n 



int 



TCk + Ti- 



(3) 



where Xij is a renormalization factor determined by the 



expectation value of Si and Sj; and n^, 
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Cii, which 



denotes the number operator. Note that U is not a 
bare on-site Coulomb interaction but a renormalized one; 
U = Ui2 — Ji2 where U12 is the interorbital Coulomb in- 
teraction and J12 is the Hund's rule coupling between the 
orbitals 1/ = 1 and 2. 

In undoped state where one electron occupies a site 
on average in the model (we call this situation 
half filling hereafter), perovskite manganites commonly 
show the A-type antiferromagnetism. This state has 
strong 2D anisotropy attributed to antifecxjomagnetic 
stacking of ferromagnetic planes {xy planes ).al3'E3 When 
we dope holes, the system gradually recovers the three 
dimensionality and finally becomes an isotropic ferro- 
magnetic rp£ta.L±hmugh structural, magnetic and MI 
transitions.Elo'QBO'EJcJ However, the perfect polariza- 
tion of spins in the 2D planes persists through these 
transitions; moreover, the 2D anisotropy remains for fi- 
nite dopingp^J-J^asL-until the magnetic or the structural 
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transition.l3130lfir£ZI In this situation, the spin polarized 
model is well justified within the 2D plane without the 
aposteriori and ambiguous factor Ay. 

On the basis of these considerations, here we consider 
the spin-polarized 2D model as a starting point to dis- 
cuss physical properties of manganites at and around half 
filling. The explicit form of it is given by 
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(5) 



Here we consider only the nearest-neighbor hopping in- 
tegrals which are explicitly given by ilj — — 3to/4, ifj — 
-io/A, ijf = ifl = -(+)^/3^o/4 in the x{y) direction.!! 



The unit of hopping integrals to is defined in Sec. 
Our assumptions to derive Eq. (|^) are examined in Sec. 
11 D in detail. 



B. Jahn-Teller coupling 

The Mn'^+ ion has an instability of a distortion of 
MnOg octahedron. This requires a strong electron- 
phonon coiipling in the model, which is the so-called JT 
couplingo There are two phonon modes wbich^trongiy 
couple with orbital degrees of freedomE3'E3 In this 
work, lattice degrees of freedom are treated in a classi- 
cal or addiabatic approximation, that is, kinetic energy 
of phonons is neglected. l3e3 Then, the Hamiltonian for 
this JT coupling is given in the form 



cl — ph 



-g E ' 



(6) 



where g is the electron- phonon coupling constant; u and 
I are two-component vectors defined as 



Ui — Ui (cos26'i, sin20i) 
7. = 2(7^^I;-), 



(7) 
(8) 



where Ui denotes a displacement of the oxygens surround- 
ing i-th Mn site and Oi is a coupling angle defined later 
in Eq. T^{^, = x,y,z) is a pseudo-spin operator 

which denotes two degrees of freedom of the degenerate 
orbitals as 



1 

2 ^ 



'^T^ c' c . 



(9) 



where r is the Pauli matrix. 

The term given by Eq. (|^) is rewritten by substituting 
the definition (H) as 



(10) 



The last term in the Hamiltonian 
tic energy of phonons as 



denotes the elas- 



(12) 



where fc is the spring constant. 

When an oxygen connecting i and j-th Mn sites is dis- 
placed by an amplitude u^, the hopping integral between 
Mn sites, , is modulated since it originally consists of 
a product of hopping integrals between Mn and O sites. 
Note that in both the phonon modes considered here, 
oxygens shift only in the direction of Mn-O-Mn bonds. 
The Mn-0 hopping is-anodificd as (1 ± lui) under 
the displacement of i n the length unit mentioned in 
the next subsection Sec. 



HC 



Therefore, we take the 
unit of Mn-Mn hopping integrals in Eq. (|5|) as 



tn = <n 1 - 4Mf 



-7/2 



(13) 



C. Model parameters 

We set to and the Mn-Mn lattice constant as the en- 
ergy unit and the length unit, respectively. Our Hamil- 
tonian contains the following parameters; the Coulomb 
interaction J7, the electron-phonon coupling constant g, 
the spring constant fc, and the parameters Ui and di 
which determine local IT .diatppi-ions. From experiments 
and band calculations JlJ^OOEil to is considered to be 
around 0.5eV and IJ is estimated as several eV. In this 
work, C//to is considered as a parameter in the following 
calculations. From our numerical results in Sec. [II, the 



realistic value of U is discussed and compared with that 
in other theoretical investigations in Sec. IV. The spring 



constant k is roughly estimated from the frequency of an 
oxygen bond stretching phc»imi,|as the order of 10 ^ 100 
eV; we use here fc = lOO.u'EaEll The electron-phonon 
coupling constant g is hard to estimate from available 
experimental data. It may, however, be the order of leV; 
we take g = 10 in this work. The values of fc and g 
are examined in Sec. [VA by quantitative comparison 
between our results and experiments. The parameters 
Ui and di are determined to minimize the total energy. 
For simplicity, they are treated in a mean-field scheme as 
described in Sec. IIIB2 in detail. 



where m^^ = "i\y"iiv is the number operator for the new 
operator 7 which is a linear combination of the operator 
c defined as 



7ii 

7i2 



cos Qi sin Qi 
— sin Bi cos Qi 



Cil 
Ci2 



(11) 



Then, Eq. ( pi)| ) shows that the JT coupling in the clas- 
sical approximation can be expressed as the local fields 
which lift the degeneracy of two orbitals diagonalized in 
the rotated basis by the coupling angle 6i as in Eq. (11). 



D. Energy scale and doping range 

In the following, the model given by Eq. (||) is consid- 
ered in two dimensions. Here, we discuss the energy scale 
and the region of the doping concentration S where the 
present model is justified. In particular, two main sim- 
plifications in the model (^) are examined; one is neglect 
of three dimensionality and the other is the full polariza- 
tion of spins. We also comment on the t2g-spin exchange 
interaction neglected in Eq. (||). 
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In undoped states, the largest energy scale in our con- 
sideration is the charge gap Ac of the order of leV. The 
formation of the gap is a common feature irrespective of 
the dimensionality. Therefore, physics of the charge gap 
formation in the energy scale of Ac may well be captured 
in the 2D model equally to the isotropic 3D system. 

When the energy scale is lowered below leV, the 
system becomes insulating and excitations are de- 
scribed through spin and orbital exchange processes be- 
cause single-particle charge fluctuations are suppressed. 
Through these processes, the JT distortion and the or- 
bital ordering take place in the energy scale of the order 
of O.leV. Although the interlayer structure of JT distor- 
tion and orbital ordering cannot be predicted, the energy 
scale of the distortion and the orbital ordering are again 
captured by the 2D model, because they take place in a 
planar configuration and the energy scale of these order- 
ings may primarily be estimated from such anisotropic 
structure. 

Although we employ a spin-polarized model, the spin 
configuration is not an important issue in this en- 
ergy scale because the spin exchange energy is too 
small as compared to the above-mentioned energy scales. 
Whether spins are ordered or not does not make any cru- 
cial difference in our discussions. 

With further lowering of the energy scale to the order 
of O.OleV, a shortcoming of the spin-polarized 2D model 
becomes clear because the 3D spin ordering takes place 
under a competition of different possible 3D structures. 
We of course are not able to discuss such process. How- 
ever, eventually, the A-type antiferromagnetic structure 
under the 2D anisotropic orbital ordering takes place, 
which again makes spin-polarized 2D model useful due 
to the strongly anisotropic structure with ferromagnetic 
alignment of spins in each plane at very low temperatures 
^ lOK. We emphasize that the above energy hierarchy in 
undoped states provides the underlying backbone struc- 
ture of this problem. Our model is useful in all the energy 
scales except those of the 3D spin-ordering process. 

We neglect the t2g-spin interaction in the starting DE 
model given by Eq. . The exchange coupling between 
nearest-neighbor t2g spina, is roughly estimated at the 
order of 0.001 or 0.01eV.El This energy scale is in the 
same order as the 3D magnetic transition temperature 
Tn. Actually, the t2g-spui interaction is considered to 
be one of the origiaSpef the A-type antiferromagnetism 
in undoped states.E3"E3 Following the above discussions 
on the energy scales, we consider that the inclusion of 
this interaction term will not affect the essence of the 
following results in this paper. 

Upon doping, a canted antiferromagnetic st|ale_ia_px- 
perimentally observed near the undoped statc.i3lij LilEj In 
this state, the perfect polarization of spins in each plane 
remains unperturbed, however, the direction of polariza- 
tion tilts to the z direction. The canting angle increases 
with the increase in the doping concentration S and fi- 
nally the isotropic ferromagnetic state appears with MI 
transition. This recovers the three dimensionality. When 



we discuss the order-disorder transitions of the JT dis- 
tortion and the orbital ordering, this 3D effect may be 
important. 

We here estimate energy scales of the 3D effect by the 
spin canting. One is given by the transition temperature 
Tn. In Lai_2;Sr2;Mn03, Tn does not strongly depend on 
the doping concentration and remains the order of O.OleV 
through the spin-canted phase. 

Another energy scale is kinetic energy in the z direction 
gained with the increase in the canting angle. We here efe 
timate this energy gain following the work by de Gennes.Q 
He analyzed the DE model (|l|) by mean-field approxima- 
tion to explain the origin of spin-canted phase observed 
by neutron scattering experiment for Lai_a;Caa;Mn03. If 
his formalism is applied straightforwardly, roughly speak- 
ing, the kinetic-energy gain of each hole increases linearly 
with the canting angle and hence linearly with the dop- 
ing concentration. Then, the total kinetic-energy gain is 
estimated as 



-5\ 



(14) 



for Q < 5 <3 5^ where (5^ is the critical concentra- 
tion for the ferromagnetic transition. In the cas©,pf 
Lai_a;Sra;Mn03, the transition occurs at (5m — 0.16Bj't3 
In terms of our model, {P'^Y is equal to —to. 

In order to discuss the 2D anisotropic structure of the 
JT distortion and the orbital ordering in the energy scales 
smaller than Tn or Ai5f in doped states, the three dimen- 
sionality should be taken into account. In other words, 
our 2D model may be useful to discuss these ordering 
processes in the energy scales larger than Tn or ASf 
even at finite hole doping. 

In the metallic state with further doping, if we assume 
that the MI transition is a continuous or a weak first- 
order one, there is a critical region of 5 where the system 
is controlled by the criticality of the transition. In the 
real materials, the insulators show the 2D anisotropy in 
both the orbital ordering and the JT distortion. The 
metallic state near the MI transition ip^y be influenced 
by the 2D anisotropy of the insulator.ES Our 2D model 
would serve to discuss the criticality of the MI transition. 



III. RESULTS 

This section presents numerical resul ts for the model 
(Q) in the ground state. First, Sec. Ill A describes numer- 
ical technique and d etails of computational conditions. 
Then in Sec. [II B, we show numerical results at half 
filling without and with the JT coupling^ Results for 
hole-doped cases are presented in Sec. [II C| . All the 



results are discussed in comparison with experiments as 
well as with other theoretical investigations in Sec. 
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A. Method and computational details 



1. Without Jahn-Teller coupling 



We investigate ground state properties of the model 
by titpi-ttrojection quantum Monte Carlo (PQMC) 
methodB^El If JT distortions are given, that is, if the 
parameters Ui and 6i are fixed, the Hamiltonian given 
by Eq. (||) is the Hubbard model with a generalized 
hopping and external fields. Therefore we can apply a 
well-established algorithm for the , OEdipa cy, single-band 
Hubbard model straightforwardly.E3E3cjE]l In two di- 
mensions, quantum fluctuations may be important even 
in insulating states. PQMC gives unbiased results on the 
effects of electron correlation and quantum fluctuations. 

For simplicity, the JT coupling given by Eq. is 
treated in a mean-field scheme. As described in Sec. 



Ill B 2 , we fix 9i following the experimental indications 
and find a uniform solution Ui — u to minimize total 
energy by PQMC. 

For the present model, PQMC has the so-called nega- 
tive sign problem even at half filling with g = Q, although 
this problem is completely, absftnt for the ordinary Hub- 
bard model at half fiUing.HLJ However, from practical 
point of view, we have managed to reach the ground-state 
convergence in several choices of parameter values. Phys- 
ical quantities presented in the following are obtained in 
the parameter region where the negative sign problem is 
not serious in reaching the ground state properties with 
needed accuracy. The details are discussed in Appendix 
A. 

In PQMC, ground states are projected out by operat- 
ing exp(— /37i) on a trial wave function. For good conver- 
gence about the projection /3, we employ several differ- 
ent trial wave functions depending on the cases; nonin- 
teracting ^ater determinant, unrestricted Hartree-Fock 
solutions and orbital-singlet wa^e functions stabilized 
by small dimerization in hoppings.Ej The convergence for 
these trial wave functions is discussed in Appendix A. 
Discreteness Ar of the slices in the projection (3 is taken 
as Ar^o = 0.05 which is small enough to obtain ground 
state properties within numerical errorbars. 

The 2D systems with No, = I? sites are investigated 
in the following. As shown in Appendix B, we note that 
this model defined on a finite-size lattice shows differ- 
ent behavior depending on the linear dimension L = 4n 
or 4n -|- 2 (n is an integer) as well as on boundary condi- 
tions due to the shapes of Fermi surfaces in the finite-size 
systems. In extrapolations to the thermodynamic limit 
N-^ —> oo, we pay attention to choose a suitable set of L 
and the boundary condition to reduce finite-size effects. 
In the following, we mainly show the results for a series 
oi L = An with the antiperiodic boundary condition and 
L = An + 2 with the periodic boundary condition. See 
Appendix B for details. 



B. Half filling 



First, the ground states at half filling are investigated 
in the absence of the JT coupling. The system becomes 
Mott insulating for finite values of the Coulomb interac- 
tion U. In the following, we show numerical results on 
the charge gap amplitude and the orbital ordering at this 
filling. 

Figure |^ shows the system-size dependence of the 
charge gap amplitude for various values of the Coulomb 
interaction U. Here, the charge gap is calculated from 
the shift of the cheraieal potential when we dope holes to 
the half-filled state.E3 The data for each value of IJ are 
well fitteid-jpy 1/L as in the ordinary single-band Hubbard 
model.E^Ej 

Figure ^ summarizes the charge gap amplitude in the 
thermodynamic limit obtained by the extrapolations in 
Fig. [|. For comparison, we. plot the QMC results for the 
ordinary Hubbard modelo as well as the results of the 
mean-field approximations for both models. The QMC 
data on the present model show remarkable reduction of 
the charge gap amplitude in the region of U considered 
here in comparison with the mean-field results as well as 
the QMC result in the ordinary Hubbard model. This 
means that quantum fluctuations may be more impor- 
tant in the present model than in the ordinary Hubbard 
model. Although further studies are necessary to discuss 
the value of U for the opening of the gap, increasing the 
Coulomb interaction U opens the charge gap and let the 
system be Mott insulating; however in order to reproduce 
the realistic values of charge gap amplitude (Ac ~ 2to) 
only with U without the JT coupling, it is necessary to 
take much larger interaction than U — b. 

The Mott insulating state induced by the Coulomb in- 
teraction shows a staggered orbital ordering. We calcu- 
late the orbital correlation function defined as 



(fe) = 7Vs-i^Tf7;^exp(ifc. 



(15) 



with the operator given by Eq. (^) [^i — x,y, z). For a 
finite value of U , T^{k) shows a peak at A; = Q = (tt, tt) 
as shown in Fig. ^, while Ty{k) and T^{k) do not have 
a conspicuous structure. 

The moment of this staggered orbital ordering is cal- 
culated by = [T'={Q)/Ns]^^^. Figure | plots the 
system-size dependence of the staggered moment M^. 
The data are well fitted by l/L, from which we obtain 
the results in the thermodynamic limit summarized in 
Fig. ||. We also plot in the figure the QMC result for the 
antiferromagnetic. spin-ordering moment in the ordinary 
Hubbard modeo and the mean-field results as in Fig. 
^. In small U region, the moment grows slower than the 
mean-field results as the charge gap amplitude in Fig. 
^, however it grows rapidly around U — A; this may be 
related to the uniaxial property of the moment in the 
present modelr-. Actually, in numerical calculations for 
small clusters J2a the moment in our model grows rapidly 
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and is saturated to its full value = 1/2 for large 
values of [/; on the other hand, the isotropic moment in 
the ordinary Hubbard model slowly approaches a reduced 
value because of stronger quantum fluctuations. 

For the following discussions, we consider here the 
pattern of the orbital ordering. Our model exhibits a 
staggered ordering of {\dx2-y2) + \d3z2 _r2)) / {\dx-2 -y^) — 
\d3z^_r2)) type indicated by the peak structure of T^{Q). 
If we define a single-site state by the linear combination 
of the two orbitals as 



m^)) 



+ \d 



(16) 



the ordering pattern is characterized by 0i — 
(— l)l^i'+'"i l7r/4. See Fig. |^. The previous mean-field 
studies have predicted slightly |diffierciiti)atterns from our 
results, (p, ~ (-l)l^?'+^ri^/3MBEra We will discuss 



this discrepancy as a 3D effect in Sec. IV A 2. 



2. With Jahn-Teller coupling 

Next, effects of the JT coupling at half filling are inves- 
tigated. In experiments, the cooperative JT distortions 
are observed with 2D anisotropy whose ordering pattern 
is characterized by 



^ (-l)in+r?l^/6 



(17) 



witb_a uniform amplitude of u ~ 0.035 in terms of Eq. 
(^Ej Following these results, we treat here lattice de- 
grees of freedom in a mean- field scheme for simplicity; 
we fix di as Eq. (|l7|), and optimize the total energy by 
PQMC to determine a uniform solution of u. 

Figure illustrates these procedures. We change the 
values of u and calculate the total energy by PQMC to 
find the optimized ground state. In the figures, the data 
at iVs = 00 are obtained by extrapolations by using the 
size scaling relation E^gLLl^ ^gi^ — 0°) + aL~^ based 
on the spin-wave theoryj£j'c2l as exemplified in Fig. ||. For 
each value of U, the total energy is considerably lowered 
by introducing JT distortions and the data are well fitted 

by 



Eg{u) ~ Eg{u = 0) - EjT + a{u - u*f 



(18) 



near the minimum, where Ejt is the stabilization energy 
of the JT distortion and u* is the optimized amplitude 
of oxygen displacement. The data on Ejt and u* are 
summarized in Fig. ^ and will be discussed below. 

The JT distortions reduce quantum fluctuations and 
stabilize the insulating states. We calculate the charge 
gap amplitude for the optimized states at u = u*. The 
JT coupling given by Eq. (^) breaks the particle-hole 
symmetry at half filling, therefore we calculate the charge 
gap for both electron doping (A+) and hole doping (A~). 
Here, the charge gap is calculated from the imaginary- 
time dependence of uniform single-particle Green's func- 
tions defined as G{rij = 0,t) = J2i, G^{rij = 0,t) with 



G,{n„T) = e{T){^o\c,Ar)cl\^o) 

-ei-r){^o\cUr)c,M, (19) 

where j'l'o) is the normalized ground state at half filling 
and 6(r) is the step function. The recently-developed 
stable algorithm is applied to abtain G(r) for the large 
values of the imaginary time \t\£3 The results of G(t) are 
shown in Fig. ^for the case of C7 = 3 as an example. The 
exponential fits to the tails of G(t > 0) (G(t < 0)) give 

shows 
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the values of the charge gap (A"). Figure 
the system-size dependence of the charge gap amplitude 
obtained from the fits in Fig. The extrapolation in 
the 1/ L dependence gives the data in the thermodynamic 
limit. The amplitudes of the charge gap are remarkably 
enhanced by the JT distortions from the values in Fig. 

Figure gT| summarizes our results on the half-filled 
states in the presence of both the Coulomb interaction 
and the JT coupling. We plot the stabilization energy 
of the JT distortion -Ejt, the amplitude of the oxy- 
gen displacement u*, and the charge gap amplitude A^ 
as a function of the Coulomb interaction U. The cor- 
responding values in experiments on LaMnOa are also 
shown for comparison: For the stabilization energy -Ejt, 
we shoWptlLe .ordering temperature of orbital and lattice 
~ 780K.ffiBEl The temperature is not exactly the same 
quantity as Ejt, however it is useful to compare these 
quantities at least in their orders of magnitu de. This 



point will be discussed later in Sec. IV A 1. For the 
oxygen displacemftat u* , we cite X-ray experiment which 
gives u* ~ 0.035.l2j For the charge gap amplitude, we 
show the experimentaJ value Ac ~ leV obtained by the 
optical conductivity.ll3 As shown in Figs. |ri|(a)-(c), these 
experimental results are consistently reproduced in our 
calculations for the value of C/ ~ 4 or 5, although slight 
adjustment on U, g or k is required for more quantita- 
tive argument. For instance, if we take k = 130 at [/ = 5 
and g = 10, as shown in Fig. ^l], we obtain more quan- 
titative agreement between the numerical data and the 
experimental values. These results will be discussed in 
Sec. [VAl in more detail. 

The JT coupling with the angle of Eq. (^^ changes 
the orbital-ordering pattern. Figure |l^ shows the angle 
of the pattern, defined in Eq. (p^), as a function of u. 
The angle is calculated from the moments of orbital 
polarizations, Af^ and M^, in the form 



, , , 1 fM'' 
\0i\ ^ - arctan — — 



(20) 



which is easily confirmed by the definition (y^. Note 
that the case of m = corresponds to the result in Sec. 
II B 1| , that is, \(f>i\ = 7r/4. For the large values of u, 
the angle \4>i\ approaches 7r/6 which equals to \9i\ in Eq. 
(0). In the case of C/ = 4, 5 = 10, and k = 100, the 
angle is around 7r/5 for u* ~ 0.0475. See also Fig. 
This result will be compared with experiments in Sec. 
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C. Hole doping 



2. With JT coupling 



1. Without JT coupling 

In the following, the results of doped states are pre- 
sented. Carrier doping leads to an MI transition. In the 
absence of the JT coupling, the metal near the Mott tran- 
sition shows some anomalous behavior in our 2D model. 

The first anomaly appears in the charge compressibil- 
ity. Figure |l3| shows the doping-concentration depen- 
dence of the chemical potential /x. Here, in order to re- 
duce the finite-size effects, /i is calculated from the energy; 
difference between adjacent closed-shell configurations.E^ 
The chemical potential /i has a jump aX 5 — Q due to the 
finite Mott gap and changes continuously in the doped re- 
gion. The monotonic decrease of ^ with increasing 5 sug- 
gests that there is no phase separation in our 2D model 
without the JT coupling. As shown in Fig. |l^(b), the 
chemical potential scales to near the MI transition at 
(5 = 0. This indicates that the compressibility k = dn/ d/j, 
diverges as 

K(xS~^ (21) 

toward the transition. 

Another anomaly is the critical enhancement of the or- 
bital short-ranged correlations. In the doped region, the 
orbital correlation is enhanced by the Coulomb interac- 
tion. Figure |l^ shows the doping dependence of the peak 
values of the orbital correlation function defined in Eq. 
(|l|) for various values of the interaction U. For finite 
values of U, the correlation shows a diverging behavior 
towards the undoped state where the orbital correlation 
is long-ranged. This behavior becomes more conspicuous 
with increasing the value of U. 

Figure 15 summarizes the data on the orbital corre- 
lations at [/ = 4 for various system sizes. The inverse 
peak values are plotted in this figure. The data are inde- 
pendent of system sizes even near S — 0, which indicates 
that the orbital correlations are short-ranged. As shown 
in this plot, we find a scaling relation given by 



T"(Q) oc S- 



(22) 



in the vicinity of the MI transition. This suggests that 
the orbital ordering at 5 = is immediately destroyed 
by hole doping and that the short-ranged correlation di- 
verges towards the transition in the doped region. 

These anomalous properties are in marked contrast 
with the usual MI transition from a metal to a band 
insulator. The scaling relations, Eqs. ( pl| ) and (22), 
are similar to those in the [Ordinary single-band Hubbard 
model i n two d imensions.!^ This similarity is discussed 
in Sec. IV Bl to characterize the MI transition in the 



present model. Our results will be examined to consider 
the charge incoherence in experiments. 



Here, the effects of the JT coupling in hole-doped re- 
gion are investigated. Phonons are treated by the mean- 
field approximation as in Sec. IIIB2. The mean- field 



treatment generally overestimates the JT distortion es- 
pecially in doped cases since it neglects various fluctua- 
tions. Actually, Fig. [l6|(a) shows that the cooperative 
JT distortion remains finite in a wide range of doping. 

However the stabilization energy of the JT distortion 
decreases rapidly as a function of the doping concentra- 
tion d. Figure |l^(b) shows 6 dependence of the stabi- 
lization energy of the JT distortion. The rough esti- 
mates of the cutoff energy are also shown in the fig- 
ure, below which some other energy scales not consid- 
ered in our model come into play as discussed in Sec. 
[I D ; one is the |iaagnetic transition temperature Tn for 
Lai-j^Sr^^MnOsJ^J and the other is t he effec t of 3D hop- 
pings AE'j^ given by Eq. (14). In Sec. IV B 2 , the collapse 
of the JT distortion will be discussed with this rapid de- 
crease of EjT . 

Figure |l7j(a') shows doping dependence of the ground 
state energy per site iJg at [/ = 4 in the presence of the 
JT coupling. The data show a convex region close to half 
filling. 

In order to see this behavior more clearly, we consider 
the quantity in each system size defined by 



(23) 



where A* is the slope oi Eg at S — in Fig. [T^(a) 
and represents the charge gap amplitude obtained as the 
chemical potential in the closed-shell adjacent to half fill- 
ing. Figure ^b) plots doping dependence of E* . The 
convex behavior is clearly seen near the half-filled state. 
It becomes more conspicuous with increasing the system 
size. 

The convex region of Eg as a function of doping indi- 
cates an instability to phase separation. The region of 
the phase separation is determined by a tangent drawn 
from the point at (5 = to the curve of Eg in Fig. |l^(a) or 
equivalently E* in Fig. p7|(b). States which have a higher 
energy than this line in the convex region are unstable 
and replaced by ones on the line with phase mixing. The 
phase mixture consists of the half- filled state and a doped 
state at which the tangent line comes in contact with the 
energy curve. 

In the present model, the phase separation takes place 
between the half-filled state with the cooperative JT dis- 
tortion and the doped metallic state, although it is dif- 
ficult to determine the critical doping concentration in 
this case oi U = A. Within the mean-field treatment, the 
latter doped state is still accompani ed by t he JT distor- 
tion as shown in Fig. |l^(a). In Sec. IVB5 , these results 



will be examined in detail, especially to understand the 
charge-ordered state in experiments. They will also be 
compared with other theoretical studies. 
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IV. DISCUSSION 

In this section, we discuss our numerical results shown 
in the previous section. For the half- filled states, the syn- 
ergetic effects between the Coulo mb int eraction and the 
JT coupling are clarified in Sec. [VAl. Our results are 
compared with experiments on perovskite Mn oxides as 
well as with other theoretical investigatio ns. Th e pattern 
of orbital ordering is examined in Sec. IV A 2. In Sec. 



IV B 



numerical results in the hole-doped states are dis- 
cussed in comparison with the recent experiments. The 
critical properties of the MI transition, esp ecially effects 



of orbital fluctuations, are examined in Sec. IV B 1 



gins of the experimentally-observed charge incoherence. 
Doping effects in the presence of the JT distortions are 
also discussed to understand experiments in lightly doped 
region. We consider the collapse of the JT distortions in 
Sec. IV B 2 and the phase separation in Sec. IV B 3. 



A. Undoped state 

1. Synergy of Coulomb interaction and JT coupling 

Consider first the case that we take account of only 
the Coulomb inte raction and neglect the JT coupling. 
As shown in Sec. 



IIIBl 



the system becomes Mott in- 
sulating with the staggered orbital ordering. The charge 
gap opens much slower than the mean-field results, which 
shows the mean-field approximation is not appropriate 
due to strong quantum fluctuations in the region of U 
considered here. We speculate that the somewhat com- 
plicated nesting property discussed in Appendix B may 
be related to this strong fluctuation. 

This strong fluctuation effect is also inferred from in 
the previous studies. Several mean-field analyses pre- 
dicted energetic d ggerierp jsy between different spin and 
orbital structures. EHlBEa~E3 Larger values of tj are re- 
quired than that in our estimates {U ~ 4—5) to reproduce 
experimental results because they have not taken into ac- 
count large contribution from the JT coupling. Note also 
that in LDA and LSDA calculations which generally tend 
to underestimate the electron correlation effects, a ferro- 
magnetierffletal would be realized if the JT distortion is 
absent .E3e3 

On the other hand, if only the JT coupling is taken 
into account, we obtain a cooperative JT distortion with 
an orbital ordering; however, as shown in Fig. the JT 
distortion is too small to reproduce the realistic oxygen 
distortions and the charge gap amplitude. This short- 
coming may correspond to the insufficient charge-gap 
amplitude pred icted in LSDA using the experimental JT 
distortions 

When the effects of both the Coulomb interaction and 
the JT coupling are considered as in Sec. IIIB2, the 
experimental results on charge, orbital, and lattice de- 
grees of freedom are quantitatively reproduced within the 



present model. Our finding of the synergetic effect by the 
electron correlation and the JT coupling suggests that the 
value of U around 4 or 5 would be sufficient to reproduce 
experimental results. Our results give a consistent pic- 
ture of the energetic hierarchy in the undoped state: (i) 
The charge gap of the order of leV opens as the conse- 
quence of the synergy between the on-site Coulomb inter- 
action and the atomic JT coupling. These two elements 
work cooperatively and nonlinearly. (ii) The correlations 
of the orbital and the lattice distortion grow in the en- 
ergy scale one order of magnitude smaller than the charge 
gap amplitude. The transition temperature of these or- 
derings is around O.leV. (iii) The two elements which 
are neglected in our model, that is, the three dimension- 
ality and spin degrees of freedom, may not be relevant 
to understand the experimental results considered here. 
The 3D antiferromagnetic state of A-type appears at the 
energy much smaller than the above two energy scales, 
namely, at the order of O.OleV. Orbital and lattice states 
may not be much affected by the magnetic ordering. 

In our analyses, the Coulomb interaction is considered 
in an unbiased way in PQMC, while the lattice distor- 
tions are treated in the mean-field scheme. Therefore, 
the effects of the interaction on the orbital correlations 
are fully taken into account up to the short-ranged and 
dynamical fiuctuations, whereas such fiuctuations of the 
JT distortions are neglected. This will lead to the follow- 
ing consequence: The charge gap amplitude, Ac, could 
be overestimated since the formation of a polaron is not 
considered. A doped hole in the insulator may cause a 
relaxation of the JT distortion around it to form a po- 
laron. Therefore, the charge gap may decrease due to the 
formation energy of the polaron. For the precise estimate 
of Ac under this effect, it is necessary to fully take into 
account of the local JT distortions. 

Another point to be discussed in our analyses is the 
comparison between the stabilization energy of the JT 
distortion iJjx and the transition temperature of the 
orbital and lattice orderings. To be precise, Sjt is 
the energy difference between the JT-distorted state and 
the undistorted one with the orbital ordering stabilized 
by the Coulomb interaction. It contains not only the 
intersite-correlation energy of the orbital and lattice or- 
derings but also the on-site JT polarization energy. On 
the other hand, the transition temperature Tc should be 
basically determined by tho-intersite correlations of the 
orbital and lattice orderingsa because the on-site JT po- 
larization may already be more or less stabilized well 
above Tc. Therefore, in our results, -Ejt may overesti- 
mate the energy scale of Tc . 

The above two points presumably require larger values 
of parameters, J7, g, or k than in the present work. We 
need further studies on the JT fiuctuation effects and also 
on a more appropriate choice of the parameter values. 
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2. Pattern of orbital ordering 

First, consider the pattern of the orbital ordering de- 
termined only by the Coulomb interaction. As described 
in Sec. Ill B 1 , our 2D model shows the pattern with (pi — 
(-1)1^"* +■'•-■ 1 7r/4 in terms of Eq. (|6|) , which is slightly dif- 
ferent from the previous mean-field results in 3D models 
as shown in Fig. ^. In order to understand this discrep- 
ancy, we discuss here effects of the three dimensionality 
neglected in our 2D model. In the z direction, the hop- 
ping integral has a strong orbital dependence; whose ex- 
plicit form is given by i}j — i}f — tjj = and iff = —to. 
This favors larger occupation of the d^22_j,2 orbital in 
3D systems than in 2D since only the d3j,2_r2 orbital has 
an overlap to the z direction. Therefore, we expect that 
the orbital-ordering pattern obtained in the present work, 
that is, {\d^2_^y2) + \d3^2^r2))/{\dx2_y2) - \d3^2_r2)) type 
may be modified in a way that a uniform component of 
d^z^-r^ increases. In other words, our result \(f>i\ = 7r/4 
gives a lower bound for the angle \(pi\ in the 3D models. 
The resultant ordering may look like \dz2_x2) /\dy2_z2) 
type. See Fig. ^. Actually, this type of orbital order- 
ing was obtained hY the-ii^ean- field approximation for 
3D electron systems .oS'EETEa 

When the JT coupling is taken into account, as shown 
in Fig. the lattice distortion modifies the orbital or- 
dering. Our result at [/ = 4 gives ~ 7r/5, which is 
very close to {d^x^-r^) /\d3y2_r2)-type. Since our calcu- 
lations neglect the three dimensionality, the angle 10^1 ~ 
7r/5 may give a lower bound as discussed above. 

Recently, the staggered orbital ordering in LaMnOa 
was experimentally observed by the resonant X-ray scat- 
tering techniquej^j although this experiment has not suc- 
ceeded in determining the angle unfortunately. In 
principle, the angle j^^l can be determined from the ip* 
tensity of the fundamental reflection in this experiment .lZI 
More detailed experiments are desired to compare our 
numerical results. 



B. Doped state 

1. Incoherence m charge dynamics 

Here, the critical properties of the MI transition in our 
2D model are analyze d in th e absence of the JT coupling. 
As mentioned in Sec. Ill C 1 , the metal shows two anoma- 



lous properties near the MI transition; the critical diver- 
gence of both the compressibility and the orbital short- 
ranged correlation. In the 2D ordinary Hubbard model, 
the same scaling relations as Eqs. ( |2l| ) and (p^)_,have 
been reported in detailed numerical investigations£3 The 
scaling theory based on the hyperscaling hypothesis de- 
scribes this MI transition with the novel universality class 
characterized hj z = 1/v — 4,e3 where z is the dynam- 
ical exponent and v is the correlation- length exponent. 
Recently, other anomalies in this MI transition have also 



been_£pund in the localization length in the iiis,ulating 
stata23 and the flat dispersion near k — (tt, 0).ll3 These 
two also support the scaling statement. The scaling rela- 
tions in our model, Eqs. ( ^l|) and (p2|), imply that the MI 
transition in the present model is in a similar universality 
class to the ordinary Hubbard model. 

This argument may give a clue to understand the 
strong charge incoherence nea£,[tiie MI transition ob- 
served in perovskite Mn oxidesES'L^I In the ordinary Hub- 
bard and t-J models, the optical conductivity has been 
numerically investigated-jbo clarify the strong incoherence 
in charge dynamics. OO In these models, the residual 
entropy and the spin short-ranged correlations are en- 
hanced towards the MI transition point, which may be 
the origin of such charge incoherence. The role of spins in 
these models is replaced with orbital degrees of freedom 
in the present model (^. This suggests that the orbital 
correlations induced by the Coulomb interaction may be 
one of the driving mechanisms of the experimentally- 
observed charge incoherence. Actually, we have also cal- 
culated the optical conductivity in the present model by 
the exact diagonalization of small clusters and obtained 
strongly incoherent properties as expected. These results 
will be reported elsewhere in detail.O 

Experimentally, the incoherent charge dynamics is ob- 
served not only in the 3D materials but also iuptiie 2D 
ones such as La2-2a;Sri+22:Mn207 near x = 0.4£3 This 
fact suggests that the essence of the charge incoherence in 
these materials is captured by the models of metals near 
the MI transition from strongly 2D anisotropic insula- 
tors. The scaling argument assuming such an anisotropy 
in energy dispepions has been discussed recently by one 
of the authors.CJ It implies that the system has a 2D 
anomalous dispersion with a flat fc-dependence caused 
by the critically-enhanced 2D orbital correlations. This 
argument provides one possible scenario to reconcile the 
puzzling situations in the 3D materials which were men- 
tioned in Sec. |; the specific- heat coefficient 7 remains 
small and is not critically enhancedEil while the Drude 
weight D and the discontinuity at the Fexmi-edge are 
strongly suppressed in the critical region.tJLj Our re- 
sults support the existence of the 2D anomalous disper- 
sion in the charge excitation near the Mott transition in 
the present model. To confirm this scenario further, it is 
necessary to investigate the fc-dependence of the disper- 
sion directly. This is left for future study. 



2. Collapse of JT distortion 

In real materials, the cooperative JT distortions per- 
sist up to finite doping concentration in the spin-ca pted. 
state, for instance, up t o (5 c± .1 in Lai-^^Sr^^MnOsJHJ'Ej 
In our results in Sec. 



IIIC2, although the JT distor- 



tion loses its stability rapidly with increasing the doping 
concentration, it remains finite in a wide range of doping. 
This presumably results from the mean-field treatment in 
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which the fluctuations of the JT distortions are neglected. 
The doped holes may form polarons, which leads to 
importance of the so-called breathing mode of phonons.EZi 
Therefore, the short-ranged JT fluctuations should be 
taken into account more seriously in the doped region 
than in the undoped state. To precisely estimate the 
doping concentration where the JT distortion collapses, 
it is necessary to fully consider the JT fluctuations. 

Since the distortion has strong 2D anisotropy, the three 
dimensionality may also be important to understand the 
collapse of the JT distortiom ^In Fig. |l^(b), the en- 
ergy scales discussed in Sec. |IIE1| are plotted for compar- 
ison; the magnetic transition temperature Tn and the 
kinetic-energy gain by the three dimensionality, AE^. 
Our model may not be appropriate to discuss phenomena 
happening in the energy scale comparable to them. In or- 
der to consider the JT distortion for J 0.1 where Ejt 
becomes comparable to A.E^ , we may need to include 3D 
effects. 

The above two points require further studies on doping 
effects with the JT distortions. The former, the fluctu- 
ation effects of the JT distortions, may be partly taken 
into account within the present model if we fully op- 
timize the parameters Ui and 9i. 



3. Phase separation 

Our results in the doped region with the JT coupling 
indicate the existence of the phase separation. This phase 
separation is a consequence of two factors; one is the 
critically-enhanced compressibility which exists even in 
the absence of the JT coupling as in Eq. and 
the other is the large stabilization energy Ejt which de- 
creases rapidly with increasing doping. Therefore, even 
if the JT fluctuation effects are taken into account, the 
phase separation may persist in the present model. 

In real systems, the additional elements such as the 
long-ranged Coulomb interaction might be crucial to re- 
place the phase separation by another instability, for in- 
stance, an ordering of the hole distribution. Recently, 



the charge ordering .has been reported in Lai_j;Srj;Mn03 
around x — 0.125Bjo From neutron-scattering data, 
this phase has been interpreted as a stacking of two al- 
ternating layers; one has little hole with a similar JT 
distortion to the undoped state, and the other contains 
excess holes and shows no apparent JT distortion. This 
phase is favorably compared to our result of the phase 
separation because it predicts the phase mixture of the 
JT-distorted undoped state and the doped one. The lat- 
ter doped phase may show no apparent JT distortion if 
we take account of the JT fluctuations as mentioned in 



Sec. [VB2 



doped ferromagnetic metal at low temperatures. Effects 
of the JT distoption on this phase separation have also 
been discussed. e3 In these cases, the phase separation 
originates from a change in the density of states between 
two states; the antiferromagnctic insulator stabilized by 
superexchange interaction has a narrow band while the 
ferromagnetic metal dominated by the DE interaction 
has a rather wide band. Note that in the 'simple' DE 
model, the phase separation takes place even in the ab- 
sence of the Coulomb interaction between the conduction 
electrons. On the other hand, in our model, the Coulomb 
interaction plays a crucial role to induce the phase sepa- 
ration as mentioned above. 



4- Remark 

Our numerical results may give a comprehensive sce- 
nario to understand the complicated properties in lightly- 
doped Mn oxides: Doped holes induce the JT fluctu- 
ations. The spin canting enhances the 3D motions of 
doped holes. These may contribute cooperatively to 
collapse the cooperative JT distortion which becomes 
unstable rapidly at a small doping concentration. On 
the other hand, the critical region of the MI transition 
may be wider than the critical concentration of the col- 
lapse. Note that the critical region of the Mott transi- 
tion spreads up to (5~0.2 — 0.3in the case of {/ = 5 as 
shown in Fig. 



13(b). Then, the doped state beyond the 



Recently, the phase separation has also been discussi 
in the original 'simple' DE model given by Eq. (|l|) .EZl"E2l 
In these studies, the system shows the phase separation 
between an undoped antiferromagnctic insulator and a 



collapse of the JT distortion shows the strong incoher- 
ence in charge dynamics. In these lightly-doped regions, 
charge degrees of freedom have an instability suggested 
by the phase separation in our results. At a commen- 
surate filling, for instance, at S — 0.125 = 1/8, doped 
holes may form a hyperlattice structure as indicated in 
experiments. 



V. SUMMARY 

In this work, roles of orbital and lattice degrees of free- 
dom in the extended double-exchange models have been 
investigated to understand remarkable properties of un- 
doped and lightly-doped perovskite Mn oxides. The mod- 
els contain the Coulomb interaction, double degeneracy 
of the conduction bands and the Jahn- Teller coupling. 
The ground state properties are studied under full polar- 
ization of spins with two-dimensional anisotropy by using 
the quantum Monte Carlo method. Fluctuation effects of 
Coulomb interaction and orbital degrees of freedom are 
investigated without biased approximations. 

In the undoped states, the importance of both the 
Coulomb interaction and the Jahn- Teller coupling is clar- 
ified. In the absence of the Jahn- Teller coupling, the sys- 
tem shows strong quantum fluctuations with remarkable 
suppression of the Mott gap amplitude. A realistic am- 
plitude of the Jahn- Teller coupling suppresses the fluc- 
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tuations and stabilizes the insulating state. Synergy be- 
tween the Coulomb interaction and the Jahn- Teller cou- 
pling remarkably enhances the charge gap amplitude and 
reproduces the realistic Jahn- Teller distortion. Our nu- 
merical results give us a consistent picture of the undopcd 
states, especially of the characteristic hierarchy of energy 
scales among charge, orbital-lattice, and spin degrees of 
freedom. The energy hierarchy consists of the charge 
gap ~ leV, orbital and lattice ordering energy ~ O.leV, 
and spin ordering and 3D coupling scales ~ O.OleV. The 
pattern of staggered orbital orderings has also been in- 
vestigated in detail. It would be desired to compare it 
with more detailed experiments in the future. 

Upon doping of carriers, the critical properties of the 
metal-insulator transition have been investigated in the 
absence of the Jahn- Teller distortion. Our numerical re- 
sults show two anomalies; the critical enhancement of the 
compressibility and the orbital short-ranged correlations 
towards the transition. These results are analyzed by 
using the scaling argument in comparison with the ordi- 
nary Hubbard model. The orbital fluctuation enhanced 
by the Coulomb interaction is discussed as one of the 
driving mechanisms of the strong incoherence in charge 
dynamics observed in experiments. 

Next, the effects of the Jahn- Teller coupling in the 
hole-doped region have been investigated. Our mean- 
field treatment of phonons suggests that the cooperative 
Jahn- Teller distortion becomes unstable rapidly with in- 
creasing hole doping. Moreover, the instability of phase 
separation appears near the undoped state. These re- 
sults are favorably compared with the collapse of the 
Jahn- Teller distortion and the charge-ordered phase ob- 
served in experiments. We propose a scenario to under- 
stand the rich physics of the lightly-doped Mn oxides; 
the order-disorder Jahn- Teller transition, the charge or- 
dering near the magnetic and metal-insulator transition, 
and the strong charge incoherence. More detailed study 
on fluctuation effects of the Jahn- Teller distortions as well 
as on the three dimensionality is left for further study. 
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APPENDIX A: 

In this Appendix, we discuss the negative sign problem 
in PQMC for the model given by Eq. (|). The ratio 
of the negative sign in QMC sampling depends on the 
system size TVs , the Coulomb interaction U^yihe trial wave 
function, and the doping concentration (5.l3'L3 

First, it should be pointed out that PQMC for the 
present model shows negative-sign samples even at the 
particle- hole symmetric point, that is, at half filling 
(/i = 0). In contrast to this, for the ordinary Hubbard 
model, the negative sign problem completely disappears 
at half fillingj£j'E2l the system is mapped to an attractive 
Hubbard model under the transformations as 

which ensure the positivity of samples in QMC. This 
mapping does not apply to the present model (j^) because 
of the difference in the diagonal hoppings, t^^ ^ t^^, and 
the off-diagonal elements i^^, t}^ . This leads to the neg- 
ative sign problem even at half filling. 

Figures [l^pO| show dependence of the expectation 
value of the sign of the MC samples, (5), at half filling. 
In the figures, (5 is the projection parameter; the ground 
state is projected out by the operation of exp(— /?7i) on a 
trial wave function. As shown in Figs. ^ and larger 
system size or larger interaction makes the negative sign 
problem more serious. 

Figure ^ indicates that it is important to choose an 
appropriate trial state in decreasing the sign problem. In 
the present work, the half-filled cases are studied by us- 
ing the unrestricted Hartree-Fock solutions and the non- 
interacting Slater determinant or the orbital-singlet wave 
functions depending on the values of I] . The unrestricted 
Hartree-Fock solutions are optimized with the value of 
the interaction [/ used in the Hartree-Fock calculations 
to give the best convergence as shown in Fig. ^ later .EJ 
The orbital-singlet wave functions are given by the sys- 
tem with a tiny dimerization introduced in the hopping 
integrals to ensure that the ground state is singlet. lj 

Figure shows b dependence of {S) at /? = 5. We 
note that (S) has the largest value at half filling. The 
negative sign problem becomes severe in the lightly doped 
region, < (5 ^ 0.5. For the doped cases, we use the 
noninteracting Slater determinant or the orbital-singlet 
wave functions as the trial state. 

The convergence of physical quantities with increas- 
ing /? strongly depends on the choice of trial wave func- 
tions. As an example. Fig. ^ shows the peak value 
of the orbital correlation function at half filling. In 
this fast convergence is obtained if we control 

the input value of [/ in the unrestricted Hartree-Fock 
calculations .Eil Generally, the choice of input U somewhat 
less than the value of the interaction in the model makes 
the fastest convergence. This may be related to the large 
deviations of the QMC results from the mean-field ones 
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due to the strong quantum fluctuations discussed in Sec. 



IIIBl 



APPENDIX B: 

Here, the noninteracting case {U = 0) in the 2D model 
is discussed. Wc show the energy dispersions of this 
model and Fermi surfaces at half filling which have a 
perfect nesting property. We also discuss the boundary 
conditions and finite-size effects in the linear dimension 
L of the system. These considerations are helpful for 
systematic analyses of numerical results in Sec. III. 

By Fourier transformation, the hopping matrix f^^ for 
Ui = in Eq. (|) is denoted as 

e-'^'ik) = eo{k)S,,' + ei(fc)x^,, + e2{k)Tf,^,, (Bl) 

where 5,^,^' is the Cronecker's delta function and is fi 
component of the Pauli matrix. Here we define 



Eo = — ^0 (cos kx + cos ky) 
El = —to/2 (cos kx + cos ky) 
£2 — V3t()/2 (cos kx — cos ky) . 



(B2) 
(B3) 
(B4) 



The 2x2 matrix e"^" (k) is easily diagonalized to derive 
energy dispersions in the form 



£±= eo± 



(B5) 



which are shown in Fig. 23 

A 
Fig 



At half filling, we obtain the Fermi surfaces shown in 
I. Note that they satisfy a perfect nesting con- 
dition, that is, k{e+ = 0) = fc(e_ = 0) -I- Q, where 
Q — (tTjTt). Therefore, we expect that the system may 
become an insulator when we switch on the interaction 
U. 

It should be stressed that the Fermi surface in Fig. |2j 
has many different ways of nesting. Besides the above- 
mentioned nesting with the wave vector Q between two 
different branches, it is possible to nest the Fermi sur- 
faces in the same branch with wave vectors (tt, a) or (a, tt) 
where a is an arbitrary number between and n. In ad- 
dition, the nesting between different branches with (tt, a) 
or (a, tt) is also possible. Therefore, our 2D model (|) 
may have a competition between these various nestings 
at half filling, which may lead to the strong quantum 



fluctuations as seen in numerical results in Sec. Ill B 1 



At half filling, finite-size systems with Ns ~ L sites 
are devided in two different series depending on a linear 
dimension L and boundary conditions: One is a group 
which has fc-points on the Fermi surfaces and the other 
does not have. In other words, the former has gapless 
but the latter has finite gap excitations in noninteract- 
ing cases. The former is a series in which L = An {n 
is an integer) with the periodic boundary conditions or 
L = An + 2 with the antiperiodic boundary conditions. 



and the latter is L = An with the antiperiodic bound- 
ary conditions or L = An + 2 with the periodic boundary 
conditions. 

These two series should be independently extrapolated 
to the thermodynamic limit Ng —> oo. Figure |2^ shows 
these extrapolations of the ground state energy. The ex- 
trapolated values in two series agree well within the er- 
rorbars. In Sec. Ill, numerical calculations and extrapo- 
lations to the thermodynamic limit are carefully done by 
considering these different series separately. 

In the doped cases, in order to reduce finite-size ef- 
fects, we consider mainly the electron fillings at which the 
closed-shell configuration is satisfied i n the n oninteract- 
ing caseH As seen in Sec. Ill CI and IIIC2, especially 
for the calculations of the chemical potential, this proce- 
dure works well and gives a smooth 6 dependence insen- 
sitive to the system sizes except for the phase-separated 
region. 
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FIG. 1. System-size dependence of the charge gap am- 
plitude at half filling. The circles, squares, and diamonds 
correspond to U = 3, 4, and 5, respectively. The data are 
fitted by l/L and the extrapolated values to L ^ oo are plot- 
ted as the filled symbols. The crosses show the data on the 
ordinary Hubbard model at (7 = 4£j 
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FIG. 2. U dependence of the charge gap amphtude in 
the thermodynamic limit. The circles and cross show the 
QMC data on the present model and on the ordinary Hub- 
bard model at U = Uf5 respectively. The gray and dotted 
curves are the mean-field results for the present model and 
the ordinary Hubbard model, respectively. 




FIG. 3. The orbital correlation function {x component) 
at i; = 4 in the system with A^s = 10 x 10 at half filhng. The 
data are presented along the bold line in the k space as shown 
in the inset. 
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FIG. 4. System-size dependence of the moment of the 
orbital polarization at half filling. The circles and squares are 
for [7 = 3 and 4, respectively. The data are fitted by 1/L and 
the extrapolated values to L ^ oo are plotted as the filled 
symbols. 




u 

FIG. 5. U dependence of the long-ranged order of the or- 
bital polarization in the thermodynamic limit. The circles and 
cross show the QMC data on the pessent model and for the 
ordinary Hubbard model at f/ = [/J_J respectively. The gray 
and dotted curves are the mean-field results for the present 
model and for the long-ranged antiferromagnetic order of the 
ordinary Hubbard model, respectively. 
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FIG. 6. Single-site orbital state defined as Eq. (|l|). The 
gray arrows at 0^ = and 7r/2 show the original basis in 
the Hamiltonian, d^2_y2 and d3^2_r2, respectively. The bold 
black arrows at (jji = ±7r/4 correspond to the pattern of the 
orbital ordering induced by the Coulomb interaction in our 
calculations. The other arrows show the basis at = ±7i"/6 
and ±7r/3. See the text for details. 
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FIG. 7. u dependence of the ground state energy at half 
filling; (a) (7 = 3, (b) [/ = 4 and (c) U = 5. We take g = W 
and k — 100. The circles, squares, triangles, crosses, and 
diamonds correspond to the data for L = 4,6, 8, 10, and oo, 
respectively. The data for L — oo are obtained by the ex- 
trapolation shown in Fig. ^. The curves show the fits by Eq. 
(p^) for the data for L — oo. 
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FIG. 8. System-size extrapolation of the ground state en- 
ergy for [/ = 4 at half filling. The circles, squares, diamonds, 
and triangles correspond to the data for u = 0.0, 0.02, 0.03, 
and 0.08, respectively. The data axe fitted by and the 
extrapolated values to L ^ oo are plotted as the filled sym- 
bols. 




FIG. 9. Imaginary-time dependence of the uniform 
Green's function for the ground state at U = 3, g = 10, 
and k = 100: (a) r > and (b) r < 0. The circles, squares, 
triangles, and crosses correspond to the data for L = 4,6, 8, 
and 10, respectively. Note that the origin for the log G(r) 
has an offset to distinguish the data for each system size in 
both figures. The lines are the fits to the exponential tails to 
calculate the charge gap amplitude. 
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FIG. 10. System-size dependence of the charge gap ampli- 
tude for (a) electron doping and (b) hole doping. The circles, 

squares, and triangles correspond to the data at i/ = 3, 4, and 
5, respectively. The lines are the fits of 1/L. 
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FIG. 11. U dependence of (a) the stabilization energy of 
JT distortions, (b) the amplitude of oxygen distortions, and 
(c) the charge gap amplitude. We take <; = 10 and k = 100. 
In (c), the open (filled) symbols denote the charge gap for 
the electron (hole) doping. The dotted lines show the exper- 
imental results of the corresponding quantities. The squares 
denotes the data at g = 10 and k = 130. See the text for 
details. 




FIG. 12. u dependence of the angle (j>i defined in Eq. 
(^. We take U = 4, g = 10, and k = 100. The circles, 
squares, triangles, and diamonds correspond to the data for 
L = 4, 6, 8, and oo, respectively. The lines connecting the 
data are guides to eye. 
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FIG. 13. Doping dependence of the chemical potential for 
U — 5. The open circles, squares, and triangles correspond 
to the data for L — 4,6, and 8, respectively. The solid circles 
are the data for L = 4 with the periodic boundary condition. 
The solid diamonds at 5 = show the charge gap amplitude 
obtained in Fig. |^. The data are fitted by 5^ in (b). 
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FIG. 14. Doping dependence of the peak values of the or- 
bital correlation function for A^s =6x6. The circles, squares, 
and triangles correspond to the data for U = 3,4, and 5, re- 
spectively. The dotted line shows the result for the nonin- 
teracting case. The lines connecting the data are guides to 
eye. 
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FIG. 15. Doping dependence of the inverse peak values 
of the orbital correlation function for U — 4. The circles, 
squares, triangles, and crosses correspond to the data for 
L = 4,6, 8, and 10, respectively. The gray line shows the 
fit by a (5-linear function for < 5 < 0.15. The dotted line 
shows the result for the noninteracting case. 
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FIG. 16. The JT distortion in the hole-doped region for 
U = 4, g = 10, and k = 100; (a) the oxygen distortion and 
(b) the stabilization energy. The circles, squares, triangles, 
and crosses correspond to the data for L = 4, 6, 8, and 10, 
respectively. The filled diamonds in (b) are the magnetic 
transition temperature in experiments. The gray curve shows 
the estimation of the kinetic-energy gain in the z direction 
given by Eq. (111). 
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FIG. 19. 13 dependence of the negative sign for A^s = 8x8 
at half filhng with 3 = 0. The circles, squares, and triangles 
correspond to the data for f/ = 2, 3, and 4, respectively. The 
orbital-singlet state is used as the trial wave function. The 
lines are guides to eye. 
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FIG. 17. Doping dependence of the ground state energy 
per site for U = 4, g = 10, and k = 100. The squares, 
triangles, and crosses correspond to the data for L = 6, 8, 
and 10, respectively. The lines are guides to eye. In (b), the 
quantity defined in Eq. is plotted to show the convex 

region clearly. See the text for details. 
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FIG. 20. /? dependence of the negative sign for U = 4 and 
g = in the system with Ai's = 10 x 10 at half fiUing. The cir- 
cles and squares correspond to the data on the orbital-singlet 
wave function and on the unrestricted Hatree-Fock solution 
with U — 2.1, respectively. The lines are guides to eye. 



FIG. 18. /3 dependence of the negative sign for U — 4 and 
g = at half filling. The circles, squares, triangles, and crosses 
correspond to the data for L = 4,6, 8, and 10, respectively. 
The orbital-singlet state is used as the trial wave function. 
The lines are guides to eye. 
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FIG. 21. Doping dependence of the negative sign at (3 = 5 
for U — 3 and g — in the system with A'^s = 6x6 for 
the choice of the orbital-singlet trial functions. The lines are 
guides to eye. 
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FIG. 22. Convergence for /3 of the orbital correlation func- 
tion at the peak point for 17 = 4 and jr = in the system with 
A^s = 6 X 6 at half filling. The circles, squares, and triangles 
correspond to the data obtained from the orbital-singlet wave 
function, the unrestricted Hartree-Fock solution with f/ = 2.4 
and 3, respectively. The horizontal line shows the converged 
value. 




FIG. 23. Energy dispersions for the noninteracting system 
defined in Eq. (B5); (a) £+ and (b) £_. 
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FIG. 24. Fermi surfaces in the noninteracting system at 
half filling; (a) for the branch e+ and (b) for £-. The bold 
solid lines show the Fermi surfaces and the shaded regions 
inside them are the occupied states. 
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FIG. 25. System-size extrapolation of the ground state 
energy for {7 = 3 and = at half filling. The circles are the 
data on the series of L = 4n with the antiperiodic boundary 
condition or L = 4n + 2 with the periodic boundary condi- 
tion. The squares arc the data on the series of L — An with 
the periodic boundary condition or L — An + 2 with the an- 
tiperiodic boundary condition. The data are fitted by 
and the extrapolated values to L — > oo are shown as the filled 
symbols. 



22 



